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Abstract 



We present an improved calculation of higher order corrections to the one- 
loop self energy of 2P states in hydrogen-like systems with small nuclear 
charge Z. The method is based on a division of the integration with respect 
to the photon energy into a high and a low energy part. The high energy 
part is calculated by an expansion of the electron propagator in powers of the 
Coulomb field. The low energy part is simplified by the application of a Foldy- 
Wouthuysen transformation. This transformation leads to a clear separation 
of the leading contribution from the relativistic corrections and removes higher 
order terms. The method is applied to the 2i-\/2 ^-^id states in atomic 

hydrogen. The results lead to new theoretical values for the Lamb shifts and 
the fine structure splitting. 

PACS numbers 12.20.Ds, 31.30Jv, 06.20 Jr 
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I. INTRODUCTION 



The evaluation of the one-loop self-energy of a bound electron is a long standing problem 
in Quantum Electrodynamics. There are mainly two approaches. The first, developed by P. 
Mohr , relies on a multidimensional numerical integral involving a partial wave expansion 
of the electron propagator in the Coulomb field. This approach is particularly useful for 
heavy hydrogen-like ions. The second approach is based on an expansion of the electron 
self-energy in powers of Z a, 



5EsE = -{Za)^mF, (1) 



where 



A4Q + y44i In 



{Zaf (^60 + Aei In [{Za)-^] + A^2 In' [{Za)-^] + o{Za)) . (2) 

The leading contribution as given by A^^i has been originally calculated by Bethe in 0. 
Many others have contributed to the evaluation of higher orders corrections, for details see 
an excellent review by Sapirstein and Yennie in [Q. A very general analytical method has 
been introduced by Erickson and Yennie in . Erickson and Yennie were able to calculate all 
the coefficients in (0) except for Aqq. The calculation of corrections of {Z aY relative order is 
a highly nontrivial task because the binding Coulomb field enters in a nonperturbative way, 
and there is no closed form expression for the Dirac-Coulomb propagator. Additionally, one- 
loop electron self-energy contributes to all orders in Z a, and the separation of the {Z aY 
relative contribution involves hundreds of terms. A very efficient scheme of the calculation 
has been introduced in 0. It was calculated there the Ago coefficient for the 15* and 25" 
states in hydrogen atom. The method was based on the division of the whole expression 
into two parts, and Eh-, by introducing an artificial parameter e which is a cutoff in the 
photon frequency. In the high energy part E^ one expands the electron propagator in powers 
of the Coulomb field and uses Feynman gauge. In the low energy part one uses Coulomb 
gauge and applies a multipole expansion. The most important ingredient of this method is 
the expansion in the parameter e after the expansion in Z a is performed (for details see the 
next section). 

The calculation presented in this paper is a further development of this original method. 
In the low energy part we use a Foldy-Wouthuysen transformation. The transformation 
clearly identifies the leading order contribution and separates out all higher order terms. 
An additional advantage is that the nonrelativistic Schrodinger- Coulomb propagator can 
be used here. A closed-form expression of this propagator is known in coordinate and in 
momentum space (for details see 0). This method is applied to the 2Pi/2 and 2P3/2 states. 
All coefficients including Aqq are obtained. We recover all the previously known results, and 
the new results for Aqq are in agreement with those obtained from the extrapolation of 
P. Mohr's data. Our results are relevant for single electron, small Z systems (for example 
atomic hydrogen and He^), which are currently investigated with very high precision. New 
theoretical values for the Lamb shift of the 2Pi/2 and 2P3/2 states and the fine structure 
summarize our calculations. 
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II. THE efw-METHOD 



The self-interaction of the electron leads to a shift of the hydrogen energy levels. This 
shift at the one-loop level is given by 

^^^^ = '''I (0^.^W(^l7^ ^_^_'^_ ^0^ 7.1^) - mm\^), (3) 

where 6m refers to the mass counter term, and it is understood that the photon propagator 
Dfj^iy has to be regularized to prevent ultraviolet divergences, is the Dirac adjoint = ip^ 7*^. 

For the cj-integration (/cq = ^)-, the lower part of the Feynman integration contour is 
bent into the "right" half plane with 3fJ(a;) > and divided into two parts, the low energy 
contour Cl and the high energy contour Ch, see Fig. The e parameter corresponds to the 
cut-off K which was introduced by H. Bethe in his original evaluation of the low energy part 
of the electromagnetic shift of energy levels P| (specifically, K = em). The two contours are 
separated along the line 3?(a;) = em, where e is some arbitrary dimensionless parameter, 
which we assume to be smaller than unity. This method of cj-integration has been described 
in detail in ||^. The two integrations lead to the high and low energy parts and E^, 
which are functions of the fine structure constant a and of the free parameter e. Their sum, 
however, 

SEsE{a) = EL{a, e) + Enia, e), (4) 

does not depend on e. The most important step is the expansion in e after the expansion in 
a. It eliminates, without actual calculations, many terms that vanish in the limit e ^ 0. To 
be more specific, in expanding E^ and E^ in e we keep only finite terms (the e°-coefficients) 
and the terms which diverge as e — 0. The divergent terms cancel out in the sum, the finite 
terms contribute to the Lamb shift. This cancelation of the divergent terms is an important 
cross-check of the calculation. One may use different gauges of the photon propagator for the 
two parts, because the gauge-dependent term vanishes in the limit e ^ 0. For convenience, 
we use the Feynman gauge for the high and the Coulomb gauge for the low energy part. 

In this work, the treatment of the low energy part is largely simplified by the introduc- 
tion of a Foldy-Wouthuysen (fw) transformation. It enables one to clearly separate out the 
leading (nonrelativistic dipole) term, which gives the a(Za)*^-contribution, from the rela- 
tivistic corrections, which give terms in a{Za)^. An additional advantage is the fact that all 
contributions to the low energy part can be evaluated using the nonrelativistic Schrodinger- 
Coulomb-Green's function, whose closed-form solution is well known 0. Terms which con- 
tribute to the Lamb shift up to a{Za)^ can be readily identified, and each of these can be 
calculated independently. In the low energy part we may expand in the photon momentum 
k. The terms which contribute to the Lamb shift in the order of a{Za)^ correspond to the 
"non-relativistic dipole" term (involving the non-relativistic propagator and wave function), 
the "non-relativistic quadrupole" term and the "relativistic dipole" term (which involves the 
relativistic corrections to the wave function and the Dirac-Coulomb propagator). The terms 
of higher order in k vanishes in the limit e — > 0. 

Calculations of the high-energy part are performed almost entirely with the computer 
algebra system Mathematica 0. Because of the presence of an infrared cut-off, one can 
expand the Dirac-Coulomb propagator in powers of the Coulomb potential. A subsequent 
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expansion of the propagator in electron momenta is also performed. This leads finally to the 
calculation of matrix elements of operators containing V and p on the P-states. Because 
P-wave functions vanish at the origin, all of the relevant matrix elements are finite up to 
the order of (Za)^. 



III. THE HIGH-ENERGY PART 

In this part we use the Feynman gauge {Dn^{k) — —Qn^/k^). and the Pauli-Villars 
regularization for the photon propagator 

111 

(5) 



A;2 A;2 k"^ - M'^' 
so that the following expression remains to be evaluated: 

We start by calculating the matrix element 

P = . ^ ^7^1^) (7) 

p— f — m — 

up to the order of (Za)^. The first step in the evaluation of P is the expansion of the matrix 

p— fc — m — 'y"V 

in powers of the binding field. We denote the denominator of the free electron propagator 
by D {D =^— }i — m). Realizing that the binding field V = —{ZaYm/p carries two powers 
of {Za) (with p = r/oBohr), we expand the matrix M up to V^, which leads in turn to four 
matrices, denoted Mj, 

with M = Mo + Ml + Ma + Ms + 0{(Zay). Defining Pi = (V'|Mi|V'), we write the element 
P as the sum 

P^Po + Pi + P2 + Ps + 0{{Zay). (9) 

This expansion corresponds to a division of the initial expression into 0-,l-,2- and 3-vertex 
parts. We then expand each of the matrices Mj into the standard 16 F matrices, which form 
a basis set of 4 x 4 matrices. 

15 2 

M, = 5^Q,^F^ where q,^ = - TV(F;3M,). (10) 

/3=0 ^ 
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The expansion coefficients are rational functions of the binding field, the electron and 
photon energy and momenta. They can therefore be expanded in powers of a, leaving none 
of the electron momentum operators in the denominator. Next, we evaluate the matrix 
elements of these operators with the relativistic (Dirac) wave function if). It is a property 
of P states, which vanish at the origin, that up to order {ZaY, all of the desired matrix 
elements are finite. 

As an example, we describe here the evaluation of the three-vertex matrix element P3 = 
(■j/^lMsl?/'). It takes on the same values for both 2P states. Expanding M3 into the 16 F- 
matrices, we find that up to order [ZaY, all expansion coefficients vanish except for the 
identity Id and 7°-matrices. The expansion coefficients are explicitly 

C3,id = 16K^ — — -4 = hj^V\ (11) 

[k^ + 2lj — w^) 

where /c = |k| and for simplicity m = 1, and 

^ .fc^-8A;2 + 12A;2cj + 16cu-6Ptu2-12cj2 + 4cj3-a;4_g ^ 

^3,7° = 2 .,2 I o = ^3,7» • (12) 

[k'^ + 2uj — w^) 

So up to order {Za)^, the two c-expansion coefficients are (except for their dependence on 
k and uj) functions of the binding field only. Thus, the matrix element P3 is given by 

P3 = &3,Id(^|V^'|V^) + 63,7°(V^l7V=^|V^). (13) 
The relevant matrix element of the wave function is 

{^\V'\^) = (^ItV^I^) = (Zaf + OiiZaY). (14) 

where the ffist equality holds only in the order of (Za)^. The above matrix elements take 
on the same values for the 2Pi/2 and 2P3/2 states because the radial parts of both 2P states 
are the same in the non-relativistic limit. 

For the other vertex parts, many more terms appear, and the matrix elements contribute 
in the lower order also. We give one example here, to be evaluated for the 1-vertex part, 

(^|7°P ■ (V p) 1^) = (Za^ ^'-^ (Zaf for 2P1/2. (15) 

and 

{ijb'p-{Vp)\^) = -^{Zaym'-^{Zaf for 2P3/2. (16) 

For a more detailed review of the calculations see f^. Having calculated P, we subtract the 
mass-counter-term before integrating with respect to k and u. The final k and u integration 
is performed in the following way. Those terms which appear to be ultraviolet divergent are 
regularized and integrated covariantly using Feynman parameter approach. The remaining 
terms are integrated with respect to k by residual integration and with respect to u by 
changing the integration variable to 
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u 



\/2muj — uj'^ + iuj 
\/2muj — ijj'^ — iuj 



(17) 



This integration procedure is described in details in 0. The final results for the high-energy- 
part are (for the definition of F see Eq. (|l])) 



1/2) 



+ (Za) 



4177 103, , , 103, , , 2 

In 2 In e 

21600 180 ^ ^ 180 ^ ^ 9e 



(18) 



and 



3/2) 



12 



+ {Za) 



^-^ln(2)-^ln(e)-l 
21600 90 ^ ^ 90 ^ ^ 9e 



(19) 



IV. THE LOW ENERGY PART 

In this part we are dealing with low energy virtual photons, therefore we treat the binding 
field non-pertubatively. Choosing the Coulomb gauge for the photon propagator, one finds 
that only the spatial elements of this propagator contribute. The a;-integration along Cl is 
performed first, which leads to the following expression for Ei, 

r d^k 1 
EL = -e' I 7^^™5^'^^(V'I«V'^-- ^ (^=|k|). (20) 

Hd denotes the Dirac-Coulomb-Hamiltonian, 6^ is the transverse delta function, and a' 
refers to the Dirac a-matrices. In the matrix element 

P^^ = (V^laV^-^- ^ -a^e-'^^-^l^) (21) 

we introduce a unitary Foldy-Wouthuysen transformation U, 

P" = m.\(U--e-'''U-) ^^^^_^'^^_^^^^^ (Ua'e-"''U-)m,). (22) 

The lower components of the Foldy-Wouthuysen transformed Dirac wave function ip vanish 
up to {Za)'^, so that we may approximate \Uip) by 

\U^) = 10) + \S(f)) with {(f)\6(f)) = 0, (23) 

where |0) is the nonrelativistic (Schrodinger-Pauli) wave function, and \6(j)) is the relativistic 
correction. 

We define an operator acting on the spinors as even if it does not mix upper and lower 
components of spinors, and we call the odd operator odd if it mixes upper and lower com- 
ponents. The Foldy-Wouthuysen Hamiltonian consists of even operators only. For the upper 
left 2x2 submatrix of this Hamiltonian, we have the result 0] 

H^^ = U{Hd-{E^-uj))U+ = m + Hs + 6H, (24) 
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where Hs refers to the Schrodinger Hamiltonian, and 5H is is the relativistic correction, 



5H 



+ 



2 m? 



5{v) + 



a 



Am? 



cr L 



(25) 



Now we turn to the calculation of the Foldy-Wouhuysen transform of the operators 
a' exp (k ■ r). The expression U a* exp (ik ■ r) is to be calculated. Assuming that = |k| 
is of the order 0{{ZaY), we may expand the expression U a^e^^''^U 
The result of the calculation is 



a 



1 +i(k • r) 



+7°^(l + z(k 
m 



-7 



2m3 



1 
2 

r) 
1 



(k 



2 m' 



;P [a ■ p 



in powers of (Za). 

(26) 



2^2 J.3 



- - (k ■ r)' 

2 ^ ' 

^ fr X S) 



2 m 



1 



fk ■ r) fk X - 



2 m 



7° (k X E) 



In the limit e — the odd operators in the above expression do not contribute to the self 
energy in {Z aY relative order, so one can neglect the odd operators. It can be shown easily 
that also the last term in the above expression (proportional to k x S) does not contribute 
to the Lamb shift in {ZaY relative order for e ^ 0. 

Because we can ignore odd operators, and because the lower components of the Foldy- 
Wouthuysen transformed wave function vanish, we keep only the upper left 2x2 submatrix 
of Eq. (P^, and we write U e^^'^U^ as 



p 
m 



1 + z (k ■ r) 



1 



(27) 



a 



2 m? 



2 m^' 

+_L (k • r) fk X ct'' 



r X cr 



2m 



This can be rewritten as 



m 



(28) 

where 6y^ is of order (Za)^. It is understood that the term 2- e*''''" is also expanded up to 
the order (Za)^. Denoting by E the Schrodinger energy {E = —{Za)^m/8 for 2P states) 
and by 6E the first relativistic correction to E, we can thus write the matrix element P^^ as 



(0 + 501 



p 



m 



ik ■ r 



Hs-iE-Lu) + 5H- 5E 



pj 



m 



ik ■ r 



+ 5y^ 



\4> + S4>). (29) 



In this expression, the leading term and the (first) relativistic corrections can be readily 
identified. Spurious lower order terms are not present in Eq. (^). By expansion of the 
denominator Hs — {E — uj) + 6H — 6E in powers of a, the whole expression can be written 
in a form which involves only the Schrodinger-Coulomb-Green's function 
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whose closed- form expression in coordinate space is given in Eq. (PB]). We now define the 
dimensionless quantity 

P = !^ ST,ij pij _ (31) 



Using the symmetry of the P-wave functions and Eq. (^), we easily see that P can be 
written, up to {ZaY, as the sum of the contributions ( P^ , ppj , ^ P| , |^). The leading 
contribution (the "non-relativistic dipole" ) is given by 

The evaluation of this matrix element is described here as an example. For the Schrodinger- 
Coulomb propagator, we use the following coordinate-space representation f^, 

G(ri, r2,E-uj) = J2 gii^u ^2, Yi^m (f i) (^2) , (33) 

with E — uj = —c? mj (2z/^). 

4tt? /2r,\' /2r.\' °° L^'+i f ^ii Warn /2r2\ 

av \av } \av ) ^ + l)2i+i (/ -M + - i^) 



where a = asohr = l/(am), and (A;)c is the Pochhammer symbol. The evaluation of eq. (|32D 
proceeds in the following steps: The angular integration is performed first. Secondly, the 
remaining integrals over r\ and r2 are evaluated using the formula (see e.g. ||10||), 

r( 7)r(n + /i + i)^_^ , 

n!r(/i + l) 

The following formula is useful for carrying out the summation with respect to k [[Til] , 



f M e- ^ .F. (-„, 1 + « i) . (35) 



E I ^' s»2F,(-n, 6; q 2) = r(A) (1 - s)-* 2Fi(a,6; c; --^). (36) 

n=0 is 

The summations lead to hypergeometric functions in the result, 

(3 - 6t - 3t2 + 12^3 + 29t^ + 122^5- 413 1^) 
"^"^^^ " 9(l-t)^ (1 + t)^ + 

9 (1-t)' (1 + t)' V ' vi + ty y ^ ^ 

where 



v/^27^ 1 
t = —^^^=^^^= = -p. (38) 

^-2m{E -u) 2 

In this expression the terms that gives divergent in e terms are separated out of the hyper- 
geometric function, so the could be easily integrated out. The other contributions to P (for 
definition fo P see eq. (|3T|)) are 
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the non-relativistic quadrupole 

1 



i „ik • r 



Hs-{E- Lo) 



(39) 



the corrections to the current a* from the Foldy-Wouthuysen transformation, 

1 



Psy = S^''^ {<p\6y' 



Hs-{E-tu) 



-ik ■ r I 



(40) 



the contribution due to the relativistic Hamiltonian, 



(41) 



the contribution due to the relativistic correction to the energy, 



5E 



1 



{E - cu) Hs-iE- iu] 



(42) 



and due to the relativistic correction to the wave function, 

1 



2 
3m 



Hs-{E- io] 



■P 



(43) 



For almost all of the matrix elements we use the coordinate-space representation of the 
Schrodinger- Coulomb propagator given in Eq. (0). There are two exceptions: For the non- 
relativististic quadrupole, we use Schwinger's momentum space representation and carry out 
the calculation in momentum space. A rather involved contribution is 



P5H = -—WG{E-cu) 
6m 



p) vra ^, , a 



\ m? 2 m? 



4 m? 



a- L 



G{E-uj)p^\<P). (44) 



where G{E — uj) = l/{Hs — {E — uj)). The form of 6H implies a natural separation of Psh 
into three terms. 



SH 



Pp4 + Ps + Pl.s- 



For Ps, 



Ps = -—WG{E-uj) 
6m 



27^ 



5{v) 



G{E~uo)p' 



(45) 



(46) 



which involves the zitterbewegungs-term (proportional to the 5-function), we use a 
coordinate-space representation of the Schrodinger-Coulomb propagator involving Whit- 
taker functions (this representation is also to be found in [^). The result for Psif) is 
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where 



F2{t) = 2Fjl,-2tA-2t 



t- 1 



Both terms P„4 and Pl 



s, 



1 
3m 
1 



Pp' = -—i<P\p'G{E-u;) 



8 m-^ 
a 



G{E-uo)p' 



4 m? 



a- L 



G{E-uo)p' 



involve two propagators G{E — oj). We use the Schrodinger equation and the identity 



[Hs-{E-uj),-^r] 
r or 



Za 



(47) 
(48) 

(49) 
(50) 

(51) 



to rewrite them to the form that contain only one propagator with modified parameters. 
Namely, to the desired order in (Za), the expression with two propagators can be replaced 
by an expression with just one propagator, in which an (Za) ^-correction is added to the 
angular momentum parameter / or to the fine structure constant a in the radial part of the 
Schrodinger-Coulomb propagator as given in Eq. (|33|) . For the Ppi and Pl.s contributions, 
many more terms appear in the calculation, and derivatives of the hypergeometric functions 
with respect to parameters have to be evaluated. The result consists of terms involving 
elementary functions and hypergeometric functions only, and other terms which involve 
slightly more complex functions. Some of the summations give rise to the Lerch transcendent 
$. Summations of the form 



oo a 



(52) 



can be evaluated with the help of Eq. (|5BD, for more details see Although we do not 
describe the calculations in detail, we stress that the summation with respect to the /c-index 
is the decisive point in the calculation. In general, a sensible use of contiguous relations 
is necessary to simplify the result of any of the summations. Symbolic procedures were 
written to accomplish this. Through the compartmentalization of the calculation achieved 
by the Foldy-Wouthuysen transformation, it has been possible to keep the average length of 
intermediate expressions below 1000 terms. 

The contribution to El due to the 6Ese is given by 



Er 



2a 



rem Jo 

Changing the integration variable to t, we have 

1 /-l. 1 - t2 



duuj P{uj). 



F = — dt , 



Pit). 



(53) 



(54) 
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The P-terms are integrated with respect to t by the following procedure. Terms which give 
a divergence for e — > are extracted from the integrand. The extraction can be achieved 
by a suitable expansion in the argument of the hypergeometric function(s) which appear in 
P{t). The extracted terms consist of elementary functions of t only, so they can be integrated 
analytically. After integration, the terms are first expanded in {Za) up to {ZaY, then in e up 
to e". The remaining part, which involves hypergeometric functions, is integrated numerically 
with respect to t by the Gaussian method. 

The t-integration leads to F-terms which we name according to the P-terms Pnd? -^i 
Fsy, FsH, FsE and Fs,f,. The Pnd-term, which is the same for both 2P-states, is given by 



nqj 



(55) 



We have recovered the first 9 digits of the Bethe logarithm with our (Gaussian) integration 
procedure (the value for the Bethe logarithm given in is lnfco(2P) = —0.0300167089(3)). 
The Pnd-term has, for e — *• 0, a divergence of +2/9(Za)^/e, which cancels the corresponding 
divergence in the high energy part. All other P-terms produce logarithmic divergences in 
{Zay ln(e) (see Table The results for the low-energy parts of the 2P-states are 



1/2) 



I \nko{2P) + {Zaf 



1 03 

-0.79565(1) + In 

^ ' 180 



Za) 



103 , 2 

H In (e) H 

180 ^ ' 9e 



and 

Fl{2P3/2) 



-- \nko{2P) + {Za) 



29 / 
-0.58452(1) + — ln((Za 



,-2 



90 



29 , , , 2 ■ 

— In (e) H 

90 ^ ^ 9e 



(56) 



■ (57) 



The divergence in 1/e and in ln(e) cancels out when the low- and high-energy-parts are 
added. The results for the P-factors (sum of low-energy-part and high-energy-part) are: 



P(2P, 



1/2) 



12 



In A;o(2P) + (Za) 



-0.99891(1) + ^ In ((Z 



a 



for the 2Pi/2-state and 



14 of 29 / 

P(2P3/2) = t: - 7T In fco(2P) + {Zaf -0.50337(1) + - In {{ZaY 



6 3 

for the 2P3/2-state. The Aqq coefficients are given by 

^o(2Pi/2) = -0.99891(1) 

and 



(58) 



(59) 



(60) 



^o(2P3/2) = -0.50337(1). 



(61) 



The last digit is the cumulated inaccuracy of the numerical integrations. The values for the 
A40 and Aqi coefficients are in agreement with known results 0. 
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These results can be compared to those obtained by P. Mohr [|T3| by extrapolation of his 
numerical data for higher Z, 



G5£(2) = -0.96(4), = -0.98(4) for 2Pi/2, 



and 



G'5ij(2) = -0.46(2), G5i=;(l) = -0.48(2) for 2P3/2, 
where the function Gse{Z) for 2P-states is defined by 

F = A40 + {Zaf [Ai In ({Za)-^) + Gse{Z) . 



(62) 



(63) 



(64) 

Because Gse{Z = 0) = Aqq, these values are clearly in very good agreement with the results 



of our analytical calculation. Using P. Mohr's numerical data [|12[, we have obtained the 
following estimates for higher order terms summarized by Gse,? 



F = A40 + (Zaf Aeo + Aqi In ((Za) 



{Za) GsE,7(^) 



(65) 



GsE,7(2Pi/2, Z = l) = 3.1(5) and Gse,7(2P3/2, Z = I) = 2.3(5). 



(66) 



One of the most important aspects of rather lengthy calculations such as those presented 
here is to avoid errors. The result has been checked in many ways. Except for checking the 
values of the terms divergent in e, it was also checked the value of each P-contribution as 
cij — > 0. It can be shown easily that the sum of all contributions to the matrix element P 
in the low-energy part must vanish in the limit a; — >■ 0. Care must be taken when checking 
the sum, because after the Foldy-Wouthuysen transformation, hidden terms are introduced 
which do not contribute to the Lamb shift, but contribute in the limit — 0. The hidden 



terms originate from the odd operators in Eq. (26). Taking into account these terms, the 
sum vanishes for both states. 



V. OTHER CONTRIBUTIONS TO THE LAMB SHIFT 

For the Lamb shift £, we use the implicit definition 

2 

E = mr [fin,j) - 1] - [/(n, j) - 1]' + £ + Fms, (67) 

where E is the energy level of the two-body-system and f{n,j) is the dimensionless Dirac 
energy, m is the electron mass, m.^ is the reduced mass of the system and mjv is the nuclear 
mass. 

For the final evaluation of the Lamb shift the following contributions are added: 

1. One- loop self energy. The coefficients are presented in this work. For the determination 
of the Lamb shift the reduced mass dependence of the terms has to be restored. 
The relevant formulae are given in ||^. For example, the Aqq have a reduced mass 
dependence of {mr/m)^. We use Eq. (|66|) to estimate the theoretical uncertainty from 
the one-loop contribution. 
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2. Vacuum polarization correction. It enters for P-states in higher order (for the formulae 
see II, p. 570). 



3. Two-loop contributions due to the anomalous magnetic moment []T5|. It is given in 
analogy to the one-loop contribution as 



5E. 



2— loop 



a 



TT 



m ■ 



{ZaY 



[^40 + 



(68) 



where the -B-coefficients are labeled in analogy to the A-coefficients for the one-loop self 
energy. The -B40 coefficient is due to the anomalous magnetic moment of the electron. 
It is given as 



B, 



Cji 



40 



2(2/ + 1) 



197 2, ^ 3^,^, 

+ 7r2ln2 + -C(3) 

72 6 2^^ ^ 



rrir 
m 



(69) 



whereC,, = 2(j-/)/(j + l/2). 
4. Two loop contributions in higher order. Recently, the logarithmic term 



62 



7^ — 111 { 

27 ^ 



rrir 



m 



(70) 



has been calculated in The Bq2 term, which is enlarged by the logarithm, probably 
dominates the contributions to the two-loop self energy in higher order. So the result 
may also be used to estimate the theoretical uncertainty of the two-loop contribution, 
coming mainly from the unknown Bqi coefficient. It is taken to be half the contribution 
from Bq2- 



5. Three-loop self energy as given by the anomalous magnetic moment 

3 (Za^4 



6E- 



3— loop 



a 

— I m 

TT 



[C40 + 



(71) 



where 



a 



40 



C. 



2(2/ + 1) 



1.17611('ll 



m 



(72) 



6. The additional reduced mass dependence of order {mr/mNY {Za)^ 0, which we will 
refer to as the {ZaY recoil contribution. 



5E, 



rec,4 



2n3 



N 



J- + 1/2 / + l/2^ 



'1-5, 



(73) 
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7. The Salpeter correction (relativistic recoil) in order (Zaf' as given in [Q. The formula 
is for P-states 



mmjsj TTn-^ \ S S l[l + 1){2 1 + 1) J 



(74) 



8. Relativistic recoil corrections in the order of {ZaY m^/miq, 



6E, 



m 



rec,6 



rriN 



(ZaY 



(75) 



The formula for P-states has been first calculated in This general form has been 
obtained by us. 



The above contributions are listed in table || for the 2P states. 



VI. RESULTS AND CONCLUSIONS 

The new theoretical values for the Lamb shifts of the 2Pi/2 and 2P3/2 states are 

C{2Pi/2) = -12835.99(8) kHz (76) 

and 

L{2Py2) = 12517.46(8) kHz. (77) 

From the values of the 2P Lamb shifts, the fine structure can be determined. It turns out 
that the limiting factor in the uncertatinty is the experimental value of the fine structure 
constant a. Using a value of (1987) 

a"^ = 137.0359895(61) (44 ppb), (78) 

the fine structure can be determined as 

E(2P3/2) - E{2Pi/2) = 10969043(1) kHz. (79) 

With the most recent and most precise value of a available ||18[ (1995), 

= 137.03599944(57) (4.2 ppb), (80) 

we obtain a value of 

^(2P3/2) - E{2Pi/2) = 10969041.52(9) (8) kHz, (81) 

where the first error originates from the uncertainty in a and the second from the uncertainty 
in the Lamb shift difference. Our result for the fine structure disagrees with that used by 
Hagley and Pipkin in |]19[ for the determination of L{2S — 2Pi/2)- Therefore their result of 
^(25* — 2P1/2) = 1057839(12) is to be modified and according to our calculation it should be 
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L{2S - 2P1/2) = 1057842(12) kHz. (82) 

Precise theoretical predictions for P-states could be used to compare two different kind 
of measurements of Lamb shifts in the hydrogen. One is the classic 2S'i/2-2Pi/2 Lamb shift 
measured by several groups pO|, ||2l]], |jT9|, and the second is the combined Lamb shift 



£(45* — 2S) — \jC{2S — IS") as measured by the Hansch group (for a review see |^). The 
experimental value of 2S Lamb shift can be extracted from E(2S'-2Pi/2) having the precise 
value for 2P1/2 Lamb shift, and can also be determined from the combined Lamb shift 
through the formula 



^(25) = - 



C{AS) - ^C{2S) + C{1S)) - (ciAS) - ^€(23) + £{13)) 



4 



exp \ o / theo 



(83) 



where the subscript exp denotes experimental, and the subscript theo denotes theoretical 
values. This theo combination has the property that terms scaling cancel out, which 
means that almost all QED effects do not contribute, and therefore the quantity can be 
precisely determined. Such a comparison of completely different experimental techniques is 
an interesting and valuable test of high precision experiments. 

The method of calculation presented in this paper could be directly applied for the evalu- 
ation of Lamb shifts and the fine structure in two electron systems, for example in helium or 
positronium. It was a purpose of this method to use only a Schrodinger-Coulomb propaga- 
tor, and relativistic effects are incorporated through the Foldy-Wouthuysen transformation. 
This method clearly separates out the lower and the higher order terms, and expresses the 
energy shift through the matrix elements of nonrelativistic operators. 
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Note added (2000): The analytic results for higher- order binding corrections to the Lamb shift 
0/2P1/2 and 2P^i2-states (in particular, the Aqq -coefficient) have recently been confirmed by 
an improved numerical calculation in the range of low nuclear charge numbers Z = 1-5. For 



details see the e-print [physics/0009090 
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TABLES 



contribution 


2-Pl/2 


2-P3/2 


Fnq 


-1.201150(1) + 49/90 In {e/{Zaf) 


-1.201150(1) +49/90 In {e/{Zaf) 


FSy 


0.791493(1) -2/9 In (e/(Za)2) 


0.531475(1) -2/9 In {e/{Zaf) 


FSH 


0.322389(1) -47/288 In {e/{Zaf) 


0.293749(1) - 35/288 In {e/{Zaf) 


FsE 


0.040095(1) + 5/96 In {e/{Zaf) 


0.008019(1) + 1/96 In {e/{Zaf) 


F54, 


-0.748478(1) + 13/36 In {€/{Zaf) 


-0.216612(1) + 1/96 In {e/iZaf) 


sum 


-0.79565(1) + 103/180 In {e/{Zaf) 


-0.58452(1) + 29/90 In {e/{Zaf) 



TABLE I. Contributions of relative order {Za)^ to the low energy part Fl for the 2Pi/2 ^'Hd 
2P3/2 states 



contribution 


2P1/2 in kHz 


2P3/2 in kHz 


one-loop self-energy 


-12846.92(2) 


12547.95(2) 


two-loop self-energy 


25.98(7) 


-12.79(7) 


three-loop self-energy 


-0.21 


0.10 


vacuum polarization 


-0.35 


-0.08 


{Za)^ recoil 


2.16 


-1.08 


{Zaf recoil 


-17.08 


-17.08 


{Za)^ recoil 


0.42 


0.42 


sum for 2P1/2 


-12835.99(8) 


12517.46(8) 



TABLE IL Contributions to the Lamb shift in kHz for the 2Pi/2 ^^^^ 2P3/2 states. Estimates of 
the contributions of uncalculated higher order terms are given in the text. Where no uncertainties 
are specified, they are negligible at the current level of precision. 
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FIGURES 
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FIG. 1. The (^-integration contour used in the calculation. Bending the Feynman contour Cp 
in the specified way leads to the high and low energy parts Ch and C^. Lines directly below and 
above the real axis denote branch cuts from the photon and electron propagator. Crosses denote 
poles originating from the discrete spectrum of the electron propagator. 
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